Near-surface  dynamics  of  a  separated  jet  in 
the  Coastal  Transition  Zone  off  Oregon 

A.  0.  Koch,  A.  L.  Kurapov,  J.  S.  Allen 

College  of  Oceanic  and  Atmospheric  Sciences, 

Oregon  State  University,  Corvallis,  OR,  97331 

akoch@coas .  or  egonst  ate .  edu 

Submitted  to  the  Journal  of  Geophysical  Research 
August  8,  2009 

Abstract 

Three-dimensional  circulation  in  the  coastal  transition  zone  (CTZ)  off  Oregon 
is  studied  using  a  3-km  resolution  model  based  on  the  Regional  Ocean  Modeling 
System  (ROMS).  The  study  period  is  spring-summer  2002,  when  extensive  observa¬ 
tions  are  available  from  the  North  Eastern  Pacific  Component  of  the  Global  Ocean 
Ecosystems  Dynamics  (GLOBEC-NEP)  project.  Our  main  focus  is  on  near-surface 
transports,  particularly  in  an  area  off  Cape  Blanco  where  an  energetic  coastal  cur¬ 
rent  is  separated  in  the  CTZ.  Comparisons  with  available  observations  (velocities 
from  mid-shelf  moorings,  surface  velocities  from  high-frequency  radars,  satellite  SST 
maps,  along-track  SSH  altimetry,  and  SeaSoar  hydrography)  show  that  the  model 
reproduces  qualitatively  correctly  the  flow  structure  and  variability  in  the  study 
area.  The  near-surface  flow  behavior  during  26  July— 21  August,  a  late-summer 
time  period  of  strong,  time-variable  southward  winds,  is  examined.  During  that 
period  the  coastal  jet  has  separated  from  the  continental  shelf  around  Cape  Blanco 
(43  °N).  The  energetic  separated  jet  continues  to  flow  southward  in  a  near-coastal 
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region  between  42.2  °N  and  43  °N.  It  subsequently  turns  around  42  °N  to  flow  west¬ 
ward  offshore  past  127  °W.  Relatively  vigorous  upwelling  and  downwelling  is  found 
concentrated  in  the  region  of  the  separated  jet.  Frontogenesis  secondary  circulation, 
nonlinear  effects  of  the  relative  vorticity  on  the  ageostrophic  Ekman  transport,  and 
submesoscale  instabilities  contribute  to  the  vertical  circulation  within  the  jet.  Ver¬ 
tical  velocities  are  found  to  reach  50  m  day^1  in  the  offshore  part  of  the  jet  and  100 
m  day-1  in  the  near-coastal  part,  where  the  jet  is  aligned  with  the  wind  direction. 


1  Introduction 

The  coastal  transition  zone  (CTZ)  is  a  region  of  open  ocean  adjacent  to  the  continental 
shelf  where  dynamics  are  affected  by  shelf  processes.  During  periods  of  summer  upwelling 
off  the  U.S.  west  coast,  narrow  filaments  of  cold  water  are  separated  from  the  shelf  to  the 
CTZ  [ Brink  and  Cowles,  1991].  For  instance,  off  Oregon,  a  distinctive  offshore  feature  in 
late  summer  is  a  coastal  jet  separated  near  Cape  Blanco  (42.8  °N)  [ Barth  et  al,  2000]. 
This  jet  may  reach  a  speed  of  0.8  m  s_1  at  the  surface  [Strut  et  al.,  1991]  and  carry 
cold  and  nutrient-rich  coastal  waters  as  far  as  200  km  offshore,  enhancing  dynamic  and 
biological  variability  in  the  CTZ.  This  prominent  feature  is  clearly  seen  in  satellite  SST 
imagery  (Figure  1). 

The  separated  coastal  jet  is  associated  with  various  dynamical  processes:  fronto¬ 
genesis,  nonlinear  jet-wind  interaction,  disturbances  and  instabilities  of  different  nature 
and  scale.  The  objective  of  the  present  study  is  to  investigate  the  near-surface  structure 
and  dynamics  of  the  coastal  upwelling  jet  separated  off  the  Oregon  coast  by  the  means  of 
numerical  simulations  and  dynamical  analysis. 

Despite  a  number  of  studies  that  discuss  possible  mechanisms  for  coastal  current 
separation  off  Cape  Blanco,  to  this  day  there  is  not  a  settled  opinion  on  the  dominant  dy¬ 
namical  processes.  Contributing  factors  may  include  interactions  with  topography  [ Caste - 
lao  and  Barth,  2007],  enhanced  wind  stress  south  of  Cape  Blanco  [ Samelson  et  al,  2002], 


interactions  of  the  coastal  current  with  the  southward  undercurrent  [ Barth  et  al,  2000], 
and  alongshore  pressure  gradients  set  up  during  periods  of  relaxation  from  upwelling  [Gan 
and  Allen ,  2002], 

Shelf  processes  are  relatively  better  studied  than  those  in  the  CTZ,  owing  in  large 
part  to  the  success  of  recent  coordinated  observational  and  modeling  programs,  e.g., 
Coastal  Ocean  Advances  in  Shelf  Transport  [COAST,  Barth  and  Wheeler,  2005  and  other 
papers  in  that  special  issue].  Advances  in  satellite  oceanography  have  influenced  progress 
in  understanding  near-surface  transports  in  the  CTZ.  However,  observations  of  the  three- 
dimensional  (3D)  structure  of  jets,  eddies  and  filaments  in  the  CTZ  have  been  limited. 
One  of  the  most  coordinated  efforts  in  that  regard  was  the  CTZ  program  of  1986-87  [see 
Brink  and  Cowles,  1991,  and  other  papers  in  that  special  issue].  Kadko  et  al.  [1991]  and 
Washburn  et  al.  [1991]  reported  significant  vertical  transport  of  upwelled  water  within 
the  jets  as  they  propagate  offshore.  Washburn  et  al.  [1991]  analyzed  phytoplankton 
data  to  track  coastal  water  that  moves  offshore  and  found  that  the  near-surface  water 
can  subduct  to  a  depth  of  100  m  within  the  jet  core.  The  associated  estimated  vertical 
speed  reached  6-10  m  d-1.  Dewey  et  al.  [1991]  explored  the  structure  and  dynamics  of 
a  coastal  upwelling  jet  separated  off  Point  Arena,  California  based  on  observations  made 
over  a  two  week  period  during  sustained  southward  wind.  They  found  asymmetry  in 
the  cross-jet  relative  vorticity  and  evidence  for  downwclling  and  upwelling  within  the  jet. 
The  estimated  maximum  vertical  velocities  reached  40  m  d-1  and  were  thought  to  be 
associated  with  the  observed  asymmetry  in  the  relative  vorticity  held. 

Numerical  modeling  studies  focused  on  California-Oregon  shelf  hows  [e.g.,  Gan  and 
Allen,  2002;  Castelao  and  Barth,  2007]  do  not  always  correctly  reproduce  the  westward 
extent  of  separation  in  the  CTZ,  potentially  due  to  a  limited  domain  size,  and  ideal¬ 
ized  boundary  conditions  and  wind  stress.  However,  jets  extending  oh  main  topographic 
coastal  features  were  obtained  in  multi-year  simulations  using  a  regional  scale  5-km  res¬ 
olution  model  of  the  U.S.  west  coast  run  with  seasonally-varying  atmospheric  forcing 


[ Marchesiello  et  al. ,  2003].  Springer  et  al.  [2009]  developed  a  3-km  resolution  model 
based  on  the  Regional  Ocean  Modeling  System  (ROMS)  with  realistic  atmospheric  forcing 
and  boundary  conditions  provided  by  the  Naval  Coastal  Ocean  Model-California  Current 
System  (NCOM-CCS)  [ Shulman  et  al,  2004],  The  3-km  ROMS  model  reproduces  coastal 
jet  separation  off  Cape  Blanco  qualitatively  correctly.  That  study  focused  more  on  shelf 
processes,  and  did  not  examine  dynamical  processes  in  the  CTZ,  in  particular,  the  jet 
structure.  In  this  paper,  we  utilize  a  similar  model  configuration  to  simulate  flows  in 
summer  2002,  when  extensive  observations  from  the  Global  Ocean  Ecosystems  Dynamics 
in  the  North-East  Pacific  (GLOBEC-NEP)  held  program  [ Strub  et  al,  2002]  are  available. 
We  use  that  data  to  evaluate  model  performance.  Then,  we  analyze  the  structure  of  the 
jet  separated  off  Cape  Blanco  with  an  emphasis  on  near-surface  behavior,  particularly 
under  strong  wind  conditions. 

2  Model 

Our  model  is  based  on  ROMS,  which  is  a  free-surface,  terrain-following,  primitive  equation 
ocean  model  widely  used  by  the  scientific  community  for  various  applications  [e.g.,  Haid- 
vogel  et  al,  2000;  Marchesiello  et  al,  2003;  Di  Lorenzo,  2003].  Algorithms  that  comprise 
a  ROMS  computational  kernel  are  described  in  detail  by  Shchepetkin  and  McWilliams 
[2003,  2005], 

The  model  computational  domain,  shown  in  Figure  1,  extends  from  40.5  °N  to 
47.5  °N  in  the  meridional  direction  and  from  the  coast,  near  124  °W,  offshore  to  129  °W. 
The  grid  has  approximately  3  km  horizontal  resolution  and  40  terrain-following  layers  in 
the  vertical  with  a  relatively  better  resolution  near  surface  and  bottom.  Bottom  topogra¬ 
phy  is  composed  by  merging  two  sets:  a  high  resolution  (12”)  NOAA-NGDC  Bathymetry 
of  the  U.S.  west  coast,  representing  features  on  the  shelf  and  continental  slope,  and  a 
lower  resolution  (5’)  ETOP05  product  [NGDC,  1988].  A  minimum  depth  of  10  m  is  set 
along  the  coastline. 


The  study  period  is  from  1  April  to  31  August  2002.  Initial  and  boundary  condi¬ 
tions  are  obtained  from  the  9  km  horizontal  resolution  NCOM-CCS  model  [ Shulman  et 
al. ,  2004]  that  spans  between  (134.5  °W,  116  °W)  and  (30  °N,  48.5  °N).  The  NCOM-CCS 
solution  was  constrained  by  assimilation  of  SSH  and  SST  using  a  nudging  approach  [Shul¬ 
man  et  al. ,  2004,  2007].  At  the  open  boundaries  of  our  model,  the  free  surface  elevation, 
barotropic  and  baroclinic  velocities,  water  temperature  and  salinity  are  provided  daily. 
Radiation  conditions  in  combination  with  a  relaxation  nudging  term  are  applied  for  baro¬ 
clinic  velocities  and  for  temperature  and  salinity  at  open  boundaries  [ Marchesiello  et  al. , 
2001].  Flather  conditions  [ Flather ,  1976]  are  utilized  for  the  normal  barotropic  velocities, 
and  Chapman  conditions  for  the  free  surface  [ Chapman ,  1985]. 

To  force  our  model,  we  calculate  surface  wind  stress  from  the  same  time-  and  space- 
dependent  surface  wind  velocity  fields  that  were  originally  used  to  force  NCOM-CCS, 
provided  from  the  Coupled  Ocean/Atmosphere  Mesoscale  Prediction  System  (COAMPS) 
[Hodur,  1997]  with  9  km  horizontal  and  daily  temporal  resolution.  Other  atmospheric  vari¬ 
ables  (air  temperature,  pressure,  relative  humidity,  precipitation,  and  solar  short-wave  ra¬ 
diation),  used  to  compute  atmospheric  fluxes  based  on  bulk  flux  parameterizations  [Fairall 
et  al,  2003],  are  obtained  from  the  National  Center  of  Environmental  Prediction  (NCEP) 
reanalysis  [Kalnay  et  al,  1996]  and  are  provided  as  monthly  fields,  describing  seasonal 
variability.  Figure  2  shows  time  series  of  the  wind  stress  at  two  mid-shelf  locations  within 
the  model  domain,  near  Newport  (44.6  °N)  and  the  Rogue  River  (42.4  °N),  projected  onto 
their  respective  major  principal  axes,  which  are  approximately  aligned  with  the  coast. 
During  spring-summer  the  wind  stress  over  the  Oregon  shelf  is  predominantly  southward 
with  rare  and  short  events  of  northward,  downwclling  favorable  winds.  The  wind  stress 
is  substantially  larger  south  of  Cape  Blanco  [ Samelson  et  al. ,  2002], 

The  effects  of  vertical  turbulence  are  calculated  using  the  Mellor  and  Yamada  [1982] 
2.5-level  turbulence  closure  scheme,  modified  by  Galperin  et  al.  [1988].  Horizontal  turbu¬ 
lence  is  parametrized  using  a  harmonic  term  with  eddy  diffusivity  and  viscosity  coefficients 


of  10  m2  s_1.  To  minimize  effects  of  reflection  at  the  boundaries,  a  sponge  layer  is  in¬ 
troduced  in  an  area  of  width  120  km  along  the  open  boundaries,  in  which  horizontal 
dissipation  is  gradually  increased  to  30  m2  s^1  toward  the  edges  of  the  domain. 

3  Model-data  comparisons 

3.1  SST 

A  comparison  of  the  model  monthly-averaged  SST  with  5-km  horizontal  resolution  monthly 
SST  composites  from  the  Geostationary  Operational  Environmental  Satellite  [GOES,  Ma- 
turi  et  al,  2008]  has  shown  that  the  model  simulates  the  seasonal  development  of  the 
surface  temperature  held  qualitatively  correctly  (Figure  1).  In  the  plots  of  the  satellite 
measurements  for  June  and  July,  patches  of  cooler  temperature  can  be  associated  with 
clouds.  The  extent  of  the  temperature  front  off  Cape  Blanco,  particularly  apparent  in 
July  and  August,  is  similar  in  the  model  and  in  the  satellite  imagery.  The  features  of 
the  temperature  front  in  the  monthly  means  are  sharper  in  the  model  plots  than  in  the 
satellite  plots.  We  hypothesize  that  small  scale  horizontal  eddy  fluxes,  unresolved  in  our 
model,  could  contribute  to  smearing  the  monthly  mean  observed  fronts,  and  to  increased 
temporal  variability  of  the  jet  position. 

3.2  ADCP  velocities 

Continuous  time  series  of  Acoustic  Doppler  Current  Profiler  (ADCP)  velocities  were  mea¬ 
sured  in  2002  at  three  mid-shelf  moorings:  NH10,  44.6  °N,  in  water  depth  H— 81  m  [ Kosro , 
2003],  Coos  Bay,  43.2  °N,  Pt=100  m  [Hickey  et  al,  2009],  and  Rogue  River,  42.4 °N,  H= 76 
m  [ Ramp  and  Bahr,  2008].  Mooring  locations  are  shown  in  Figure  1.  Figure  2  shows  time 
series  of  40-hour  low-pass  filtered  depth-averaged  velocities  at  the  mooring  locations, 
model  and  observations,  projected  onto  their  respective  major  principal  axes  that  deviate 
slightly  from  the  meridional  direction.  At  each  mooring,  the  model-data  correlation  co- 


efficient  (CC)  is  high  (>0.68),  and  the  root  mean  square  error  (RMSE)  is  low  (<0.14  m 
s^1).  At  the  Rogue  River  mooring,  south  of  Cape  Blanco,  the  observations  show  increased 
variability  on  the  temporal  scale  of  several  days  not  well  described  by  the  model  (see  Fig¬ 
ure  2e).  This  variability  correlates  with  that  in  the  wind  stress  (Figure  2b).  The  reason 
for  this  strong  response  to  the  wind  in  this  frequency  band  is  not  entirely  understood. 
However,  we  speculate  that  part  of  that  response  may  be  due  to  remote  forcing  south  of 
our  domain  that  is  not  totally  represented  by  the  southern  boundary  conditions.  Figure  3 
shows  time-averaged  means  and  variance  ellipses  for  the  depth- averaged  currents  at  the 
mooring  locations.  Both  the  data  and  model  reveal  larger  current  variability  in  the  along¬ 
shore  direction.  The  magnitude  and  direction  of  the  mean  current  and  the  variance  are 
very  similar  in  data  and  model.  Curiously,  the  mean  current  at  the  Rogue  River  site  inside 
the  separation  zone  south  of  Cape  Blanco  is  close  to  0,  despite  large  and  predominantly 
southward  wind  stress  at  that  location. 

Time-averaged  statistics  of  the  measured  currents  from  the  three  moorings  are  com¬ 
pared  with  corresponding  model  currents  as  a  function  of  depth  in  Figure  4.  Fairly  good 
agreement  of  the  observed  and  modeled  mean  values  and  standard  deviations  for  both 
the  larger  alongshore  v  and  smaller  cross-shore  u  components  is  found  at  NH10  and  Coos 
Bay.  At  Rogue  River,  the  signs  of  the  observed  and  modeled  mean  v  differ  below  25 
m,  but  both  are  relatively  small.  The  standard  deviations  at  Rogue  River  are  similar  in 
magnitude  to  those  at  NH10  and  Coos  Bay  with  observed  values  slightly  larger.  At  all 
of  the  moorings,  the  correlations  of  the  observed  and  modeled  along-shore  velocity  v  are 
reasonably  high  while  those  of  the  cross-shore  velocity  u  are  considerably  lower.  Likewise, 
the  normalized  root  mean  square  errors  (NR.MSE)  are  substantially  lower  for  v  than  for 
u.  The  greater  success  in  modeling  the  fluctuations  in  the  larger  along-shore  currents  v 
compared  with  the  smaller  cross-shore  currents  u  is  consistent  with  previous  model  results 
for  mid-shelf  currents  off  Oregon  [e.g.  Springer  et  al. ,  2009]. 


3.3  HF-radar  currents 


To  assess  the  accuracy  of  model  surface  currents  in  the  area  around  Cape  Blanco,  we  com¬ 
pare  the  model  and  maps  derived  from  long-range  HF  radars  [ Paduan  et  a/.,  2004;  Kosro, 
2005].  HF  radar  locations  are  shown  in  Figure  1.  The  radial  velocity  components  with  6 
km  along-beam  resolution  were  low-pass  filtered,  recomputed  into  zonal  and  meridional 
components,  daily  averaged  and  mapped  on  a  6  km  regular  grid  by  P.M.  Kosro  (OSU). 
Figure  5  shows  observed  and  model  monthly  averaged  surface  current,  speed,  and  the 
RMS  speed  of  deviations  from  the  mean,  defined  as 


Urms  — 


( u  —  u )2  +  (v  —  vy 


(1) 


where  the  over-bar  denotes  a  time-average. 

The  model  and  observed  current  patterns  are  qualitatively  comparable.  In  May 
the  surface  jet  is  already  separated  from  the  Cape  and  flows  southward.  By  July,  the 
observed  currents  form  two  jets,  one  south  and  one  north  of  Cape  Blanco.  A  similar 
structure  is  seen  in  the  model,  although  it  appears  earlier  in  June.  The  observed  jets  are 
wider  and  less  energetic  than  modeled.  This  may  come  partly  from  the  smoothing  effect 
of  the  mapping  procedure  performed  on  the  radial  data  and/or  from  small  scale  eddy 
variability  unresolved  in  the  model.  By  August,  when  the  model  jet  turns  westward  at 
42  °N,  the  data  show  a  similar  westward  orientation  of  the  flow,  although  the  observed 
monthly  mean  is  diffused  over  a  larger  area.  The  magnitude  of  the  monthly  averaged 
model  current  within  the  jets  is  similar  to  observed  magnitudes  in  May-June  (0.3-0. 5  m 
s^1),  and  is  30-40  %  higher  in  the  model  (0.6-0. 7  m  s-1)  than  in  the  observations  (0.2-0. 5 
m  s^1)  in  July- August.  The  variability  of  the  current,  showed  by  Urms  (Figure  5),  is 
of  the  same  order  (0.1-0.25  m  s-1)  in  the  observations  and  the  model  in  May-July,  with 
increased  variability  found  over  the  shelf  and  in  the  CTZ  jets.  Although  variability  in  the 
HF-radar  data  in  August  (0.1-0. 2  m  s^1)  is  lower  than  in  the  model  (0.02-0.3  m  s-1),  it 


is  distributed  more  evenly  over  the  area.  This,  together  with  the  fact  that  the  observed 
mean  currents  are  more  spatially  uniform  than  the  modeled  currents,  suggests  that  the 
position  of  the  observed  separated  jet  fluctuates  over  the  domain  more  than  that  of  the 
modeled  jet.  Again,  we  speculate  that  eddy  variability  on  scales  <10  km,  not  represented 
in  our  model,  contributes  to  jet  instability  and  variability. 

Despite  the  fact  that  the  jet  is  not  seen  clearly  in  monthly  averaged  observed  fields 
in  August,  it  may  be  readily  identified  in  daily  plots  (Figure  6).  The  dates  in  Figure  6 
are  chosen  so  that  satellite  SSH  observations  are  available  along  the  track  that  passes 
through  the  area  (shown  as  the  white  line,  see  discussion  in  Section  3.4).  The  model  and 
observed  jets  in  the  daily  plots  (Figure  6)  are  found  to  have  comparable  speeds  (0.6-0. 8 
nr  s'1),  and  across-jet  spatial  scales  (10-30  km). 

3.4  SSH 

In  summer  2002,  the  AVISO  (Analysis,  Validation  and  Investigation  of  Satellite  Oceanog¬ 
raphy,  www.aviso.oceanobs.com)  satellite  SSH  altimetry  (sea  level  anomaly)  is  available 
along  the  orbits  of  the  TOPEX/Poseidon  satellite  with  10  day  periodicity  [Fu  et  al,  1994], 
Six  tracks  cross  our  model  domain  area.  Here  we  show  comparison  with  the  data  from 
track  206,  which  passes  through  the  area  covered  by  the  HF  radar  (see  Figure  6).  For  this 
analysis,  along-track  means  are  taken  out  from  both  the  model  and  observational  lines. 
Variability  in  SSH  is  of  similar  magnitude  and  horizontal  scale  in  both  the  satellite  data 
and  the  model  and  can  be  associated  with  jets  and  eddies  in  the  CTZ.  The  SSH  gradient 
is  proportional  to  the  surface  geostrophic  current  normal  to  the  track.  Using  the  satellite 
data,  the  estimated  geostrophic  current  in  the  jets  can  be  as  large  as  0.6-0. 8  m  s_1.  How¬ 
ever,  the  location  and  intensity  of  individual  jets  and  eddies  in  the  along-track  data  and 
model  do  not  necessarily  coincide.  It  is  possible  that  assimilation  of  SSH  data  in  a  model 
of  this  class  can  improve  the  representation  of  the  time-dependent  eddy  dominated  flows 


in  the  CTZ. 


3.5  Density 


In  the  beginning  of  August  2002,  vertical  cross-shore  sections  of  potential  density  erg 
were  measured  as  part  of  a  SeaSoar  survey  [O’Malley  et  al.,  2002;  Barth  et  ai,  2005]. 
The  comparisons  of  three  cross-shore  density  sections  measured  along  41.9  °N,  43.5  °N, 
44.25  °N  in  the  beginning  of  August  with  sections  sampled  at  the  same  times  and  locations 
from  the  model  results  are  presented  in  Figure  7.  There  is  a  qualitative  and  quantitative 
agreement  between  data  and  model  in  the  two  sections  north  of  Cape  Blanco  (Figure  7  a- 
d),  especially  in  the  shape,  slope  and  spacing  of  the  isopycnals  that  provide  information  on 
properties  of  the  upwclling  density  front,  and  the  associated  normal  geostrophic  flow.  The 
model-data  agreement  in  the  southernmost  section  is  not  as  close.  In  particular,  the  eddy 
shown  by  the  uplift  of  isopycnals  near  125.3  °W  (Figure  7e)  is  only  weakly  represented  in 
the  model  section  (Figure  7f).  This  behavior  reflects  the  difficulty,  found  also  in  Section 
3.4  in  connection  with  model/data  SSH  comparisons,  of  deterministically  modeling  the 
individual  filaments  and  eddies  in  the  energetic  separated  flow  region  offshore  south  of 
Cape  Blanco. 

4  Dynamical  analyses  of  jet  structure 

4.1  Mean  near-surface  circulation 

In  order  to  explore  the  mean  near-surface  circulation  in  the  CTZ  off  Oregon,  we  choose  a 
late  summer  27  day  time  interval  from  26  July  to  21  August  when  the  wind  over  the  study 
area  is  generally  upwclling-favorablc,  but  varies  in  time  on  a  typical  several-day  time  scale 
(Figure  8)  and  the  separated  jet  is  well  developed  and  extends  offshore  as  far  as  200  km 
(Figure  9).  The  region  of  relatively  large  horizontal  gradient  of  the  time-mean  surface 
density  held,  along  with  the  SSH  held,  shows  the  upwelling  front  location  (Figure  9, 
left).  The  coastal  upwelling  jet  is  rather  discontinuous  along  the  front  and  breaks  into 
a  few  separation  zones.  The  separated  jet  intensifies  and  reaches  maximum  values  of 


0.6-0. 7  m  s-1  south  of  Cape  Blanco  (Figure  9,  right).  The  maximum  jet  variability  is 
found  within  these  separation  zones.  North  of  the  Heceta  Bank  complex,  around  45  °N, 
the  separated  jet  has  approximately  equal  strengths  in  its  along-  and  normal  to  coast 
components.  In  contrast,  offshore  in  the  separation  zone  south  of  Cape  Blanco  the  jet  is 
noticeably  stronger  than  its  upstream  along-coast  link,  whose  strength  may  be  weakened 
by  the  cyclonic  mesoscale  eddy,  located  west  of  the  jet  and  centered  near  43  °N,  125.7  °W, 
that  is  stationary  through  the  averaging  interval. 

Horizontal  fields  of  the  mean  and  standard  deviation  of  the  vertical  velocity  at  25 
nr  depth  over  this  time  period  (Figure  10)  show  interesting  behavior.  In  particular,  the 
spatial  pattern  of  the  standard  deviations  of  w,  which  can  reach  magnitudes  typically  two 
times  greater  (around  20  m  day-1)  than  the  magnitudes  of  the  mean  values  (around  10  m 
day-1)  show  relatively  large  values  strongly  concentrated  in  the  region  of  the  separated 
jet.  That  region  extends  along  the  jet  from  a  location  within  about  60  km  of  the  coast 
near  Cape  Blanco  (43  °N),  where  the  separating  jet  is  directed  southward,  offshore  to 
about  127  °W,  where  the  jet  is  flowing  westward.  The  magnitude  of  the  mean  values  of  w 
(Figure  10)  are  appreciable  in  a  similar  spatial  region,  but  with  greater  values  in  the  near- 
coastal  location  between  42.2  °N  and  43  °N,  where  the  jet  flows  southward.  The  behavior 
in  that  location  is  characterized  by  concentrated  mean  upwelling  velocities  offshore  in  the 
jet,  with  downwelling  velocities  immediately  inshore. 

The  horizontal  spatial  structure  of  the  mean  and  standard  deviation  of  the  vertical 
velocity  at  25  m  depth  in  Figure  10  is  remarkable.  In  particular,  it  reflects  the  presence  of 
energetic  near-surface  ageostrophic  processes  that  are  relatively  localized  in  the  vicinity  of 
the  separated  coastal  jet  off  Cape  Blanco.  It  also  provides  motivation  for  the  analyses  that 
follows  of  the  time-dependent  dynamics  that  leads  to  the  large  variability  in  near-surface 
vertical  velocities  in  this  region. 


4.2  Lagrangian  analysis  of  surface  flows 


To  get  an  idea  of  the  general  character  of  Lagrangian  flows  in  the  CTZ  during  the  cho¬ 
sen  late  summer  time  period,  we  release  65  model  surface  particles  simultaneously  on  1 
August-0000  (UTC,  hereinafter)  along  the  81  m  isobath  (the  depth  of  mooring  NH10) 
every  10  km  in  the  meridional  direction  between  latitudes  41  °N  and  47  °N.  Model  tra¬ 
jectories  are  obtained  by  integration  of  surface  velocities,  saved  every  4  hours,  using  a 
4th  order  Runge-Kutta  method.  The  particles,  advected  by  the  surface  current,  begin  to 
group  in  offshore  directed  filaments  over  the  initial  period  1-10  August  (Figure  11).  The 
speed  of  the  particles  entrained  in  filaments  can  be  several  times  higher  than  the  speed 
associated  with  the  offshore  Ekman  transport.  For  the  trajectories  over  1-20  August 
(Figure  11)  we  can  see  two  major  compact  jets  that  originate,  respectively,  over  the  shelf 
near  the  Heceta  Bank  complex  (44  °N)  and  near  Cape  Blanco  (43  °N).  The  pattern  of 
surface  trajectories  over  the  Oregon  CTZ  (Figure  11)  is  qualitatively  consistent  with  the 
shape  of  the  observed  SST  front  in  August  (Figure  1). 


4.3  Surface  strain  rate  field 


Although  Lagrangian  tracking  described  above  is  a  powerful  visualization  tool,  it  is  lim¬ 
ited  to  areas  where  particles  are  seeded  at  an  initial  time  or  found  later.  To  quantify  the 
rate  of  relative  surface  particle  separation  over  the  entire  domain,  the  surface  strain  rate 
tensor  can  be  computed: 
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where  (x,y)  are  local  Cartesian  coordinates  [e.g.,  Batchelor ,  1967,  Section  2.3].  The  di- 


agonal  elements  of  the  matrix,  exx  and  syy,  represent  the  normal  strain  rate  and  the 
non-diagonal  elements,  eyx  =  exy ,  the  shear  strain  rate.  At  each  location,  the  strain  rate 
tensor  can  be  rotated  to  principal  axes,  in  which  the  shear  component  is  zero  and  the 
normal  components  eq  and  £2  are: 
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The  first  term  on  the  right  hand  side  of  (3)  is  the  divergent  component.  It  provides  an 
estimate  of  the  relative  particle  separation  rate  due  to  the  relative  change  of  the  surface 
area.  The  second,  deformation  component  quantifies  the  particle  separation  rate  due  to 
the  change  in  the  shape  of  a  small  surface  domain,  without  changes  of  its  area. 

Figure  12  shows  the  two  terms  in  the  principal  surface  strain  rate  held  (3),  averaged 
over  the  27  day  period,  26  July-21  August.  In  these  averaged  maps  the  deformation 
part  of  the  principal  strain  rate  is  dominant.  Note  that  there  is  divergence  in  the  surface 
current  near  the  coast  associated  with  surface  Ekman  transport,  although  in  the  monthly 
averaged  plot  its  magnitude  is  smaller  than  that  of  the  deformation  part  of  the  strain 
rate.  In  the  following  we  will  examine  the  how  during  the  wind  event  around  1  August. 
On  this  date,  offshore  between  125  °W  and  127 °W,  an  energetic  separated  jet  is  directed 
westward  approximately  along  42  °N  (Figure  13,  right).  The  wind  stress  in  the  area  of  the 
jet  is  strong  (-0.35  N  m-2)  (Figure  8)  and  is  directed  southward,  nearly  perpendicular  to 
the  jet  direction  (Figure  13,  left).  Velocities  within  the  jet  reach  0.9  m  s-1.  During  this 
event  the  divergence  term  of  the  principal  strain  rate  becomes  comparable  in  magnitude  to 
the  deformation  term,  particularly  along  the  paths  of  the  CTZ  jets  (Figure  12c, d).  Both 
terms  are  increased  in  magnitude  compared  to  the  27  day  averages  (note  that  color  bar 
limits  for  Figure  12a, b  and  Figure  12c, d  are  different).  Since  the  divergence  of  the  surface 
horizontal  how  is  entirely  due  to  its  ageostrophic  component,  i.e.  V/,-u  =  dua/dx+dva/dy 


(Section  4.4),  this  analysis  provides  an  additional  indication  of  the  important  role  of  time- 
dependent  ageostrophic  processes  in  the  near-surface  dynamics  of  the  separated  jet  south 
of  Cape  Blanco. 

4.4  Structure  of  the  separated  jet  on  1  August 

To  study  further  the  nature  of  the  flow  that  leads  to  the  large  mean  and  RMS  near¬ 
surface  vertical  velocities  in  the  region  of  the  separated  jet  (Figure  10)  we  examine  some 
instantaneous  horizontal  fields  and  vertical  sections  on  1  August,  during  a  period  of  strong 
southward  winds  (Figure  8)  when  the  offshore  separated  jet  is  well  developed  (Figure  13, 
right).  The  wind  stress  field  on  1  August  (Figure  13,  left)  shows  the  known  spatial  increase 
in  magnitude  of  the  wind  stress  south  of  Cape  Blanco  [Samelson  et  al.,  2002],  It  also  shows 
that  on  1  August  the  stress  vectors  have  a  predominantly  north-south  direction  in  the 
region  of  the  separated  jet  (41.5  —  43  °N,  124.5  —  127  °W).  We  look  first  at  the  vertical  jet 
structure  along  two  sections,  one  near-coastal,  east-west  section  along  42.63  °N  and  one 
offshore,  north-south  section  along  126  °W  (locations  shown  in  Figure  13).  Instantaneous 
values  of  potential  density  ae,  the  respective  geostrophic  along-jet  velocity  components 
vg  =  (' l/fp0)(dp/dx )  and  ug  =  —(l/fpQ)(dp/dy),  the  ageostrophic  va  =  v-vg,  ua  =  u-ug 
velocity  components,  where  /  is  the  Coriolis  parameter,  p0  is  the  reference  density  and  p 
is  the  pressure,  are  shown  in  Figure  14. 

At  the  location  of  the  near-coastal,  east-west  section,  the  coastal  jet  has  separated 
offshore  of  the  continental  shelf,  but  is  still  flowing  southward  in  the  direction  of  the  wind 
stress  (Figure  13).  The  dominant  along-jet  v  velocities  are  in  geostrophic  balance  with 
a  density  field  characterized  by  a  strong  upwelling  frontal  structure  (Figure  14).  There 
is  some  augmentation  of  the  southward  geostrophic  vg  velocities  by  an  ageostrophic  va 
component  in  the  jet  core.  The  across-jet  and  vertical  ageostrophic  velocities  exhibit 
characteristics  of  frontogenesis  secondary  circulation  (FSC)  [ Hoskins ,  1982;  Capet  et  al. , 
2008b] .  That  circulation  is  dominated  by  vigorous  downwelling  on  the  inshore,  high  den- 


sity  side  of  the  jet,  concentrated  in  a  region  with  small  horizontal  scale  of  about  10  km. 
Correspondingly  vigorous  upwelling  occurs  adjacent  on  the  lower  density  side  over  a  sim¬ 
ilar  short  horizontal  scale.  The  vertical  scales  of  the  larger  downwelling  and  upwelling 
velocities  extend  to  about  90  m  depth.  The  near-surface  ageostrophic  ua  velocity  com¬ 
ponent  is  primarily  negative  and  directed  offshore,  presumably  driven  by  the  southward 
wind  stress  in  general  accordance  with  Ekman  dynamics.  In  the  vicinity  of  the  front,  how¬ 
ever,  ua  weakens  considerably  and  reverses  sign,  as  it  takes  part  in  the  vertical  circulation 
processes.  The  basic  characteristics  of  this  wind-intensified  FSC  are  in  agreement  with 
that  predicted  by  recent  theoretical  and  modeling  studies  [e.g.,  Thomas  and  Lee,  2005]. 
We  note  that  ua  appears  to  be  larger  in  magnitude  offshore  on  the  negative  vorticity  side 
of  the  jet  which  would  be  consistent  with  the  nonlinear  effect  of  the  geostrophic  relative 
vorticity  dvg/dx  on  the  Ekman  dynamics  [Stern,  1965,  Niiler,  1969,  Thomas  and  Lee, 
2005],  represented  as 
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At  the  location  of  the  offshore  north-south  section,  the  separated  jet  is  flowing  west¬ 
ward  (Figure  13).  During  this  strong  wind  event  (and  during  the  entire  averaging  period 
26  July-21  August)  the  separated  jet  advects  cold  dense  water,  upwelled  near  the  coast, 
westward  (Figure  13).  By  the  thermal  wind  balance,  the  vertical  shear  in  ug  (Figure  14) 
is  of  opposite  signs  on  the  two  sides  of  the  jet.  As  a  result,  the  zonal  along-jet  geostrophic 
velocity  ug  has  an  asymmetric  structure  in  the  across-jet  north-south  section.  The  neg¬ 
ative  along-jet  ug  velocities  are  strengthened  in  this  section  by  a  negative  ageostrophic 
ua  component,  which  appears  to  be  primarily  wind-driven,  but  is  again  relatively  large 
locally  in  the  jet  core.  The  ageostrophic  va  component  is  negative  (southward)  near  the 
surface  and  positive  at  depths  greater  than  about  20  m,with  notably  larger  magnitudes  on 


the  northern,  negative  vorticity,  side  of  the  jet.  Energetic  vertical  circulation  is  present, 
involving  upwelling  concentrated  in  a  small  0(10  km)  horizontal  scale  region  on  the  north¬ 
ern,  lower  density  side  of  the  jet  with  adjacent  downwelling  on  a  similar  horizontal  scale 
on  the  southern,  higher  density  side.  This  vertical  circulation  appears  to  be  primarily 
associated  with  along-jet  submesoscale  instabilities  (to  be  discussed  in  Section  4.5).  We 
note  that  the  structure  of  the  vertical  velocity  is  qualitatively  similar  to  the  conceptual 
picture  of  vertical  processes  in  an  observed  CTZ  jet  based  on  the  analysis  of  hydrographic 
data  in  the  CTZ  held  experiment  [ Dewey  et  al,  1991]. 

To  examine  other  dynamical  features  associated  with  the  separated  jet  in  these 
two  locations  we  plot  in  Figure  15  corresponding  across-jet  vertical  sections  of  buoyancy 
frequency  N2  =  — ( g/po)(dag/dz ),  turbulent  vertical  diffusivity  coefficient  Kh,  turbulent 
kinetic  energy  (TKE),  including  contours  of  the  Richardson  number  Ri  =  N2  /  S2,  where 
S2  =  (du/dz)2  +  (dv/dz)2  is  the  sum  of  the  squared  vertical  shear  of  the  horizontal  velocity 
components,  and  turbulent  shear  production  P  =  KmS2,  where  Km  is  the  turbulent 
vertical  viscosity  coefficient  [e.g.  Wijesekera  et  al,  2003].  In  both  sections,  large  values 
of  buoyancy  frequency  (Figure  15a, e)  reflect  stable  vertical  stratification  and  correspond 
to  the  areas  of  relatively  large  vertical  density  gradients  (Figure  14a,f).  Areas  of  unstable 
stratification  with  negative  N2,  marked  by  the  white  stars,  result,  through  the  turbulent 
closure  scheme,  in  accompanying  large  values  of  vertical  diffusivity  Kh  (Figure  15b, f).  In 
the  coastal  section,  the  unstable  region  with  large  Kh  at  about  124.9  °W  appears  to  be 
related  to  the  downwelling  circulation  at  that  location,  but  unstable  conditions  are  found 
at  other  locations  away  from  the  frontal  regions  of  strong  vertical  circulation  as  well.  In 
both  sections,  regions  of  relatively  large  TKE  are  found  in  the  surface  layer  on  the  light 
side  of  the  front  (Figure  15c, g),  reflecting  the  results  of  relatively  large  shear  production 
P  (Figure  15d,h)  in  those  locations.  In  the  offshore  section,  the  region  of  large  surface 
layer  TKE  is  considerably  greater  and  extends  farther  from  the  jet  core.  An  additional 
analysis  (not  shown  here)  shows  that  the  enhanced  shear  production  in  that  region  is 


caused  by  the  shear  in  the  ageostrophic  across-jet  velocity  component  va. 

The  sections  of  TKE  and  K h  (Figure  15)  show  significant  spatial  variability  that 
evidently  has  some  relation  with  the  presence  of  the  separated  jet.  To  examine  that  be¬ 
havior  further,  we  plot  in  Figure  16  horizontal  fields  on  1  August  of  relevant  variables, 
including  ageostrophic  surface  velocity  vectors  superposed  on  the  relative  vorticity  ( g  of 
the  geostrophic  surface  currents,  the  maximum  value  in  the  upper  25  m  of  TKE,  and 
the  maximum  value  in  the  upper  25  nr  of  the  turbulent  closure  scheme  stability  function 
Gh  =  min(— l2 N2 /2TKE,  0.028),  where  l  is  a  turbulent  length  scale  [e.g.  Wijesekera 
et  al,  2003].  A  relationship  between  the  surface  ageostrophic  currents  and  the  surface 
geostrophic  vorticity  is  apparent  in  the  horizontal  holds.  The  direction  of  the  surface 
ageostrophic  velocity  vectors  is  locally  changed  in  areas  of  jet  hows  where  the  vorticity 
changes  its  sign  (Figure  16,  left).  The  across-jet  sections  (Figure  15c, g)  indicate  that  the 
TKE  is  strongly  affected  by  the  ageostrophic  processes  in  the  surface  layer.  The  horizontal 
held  of  the  maximum  value  of  TKE  in  the  top  25  m  (Figure  16,  middle)  gives  an  assess¬ 
ment  of  the  spatial  extent  of  that  effect.  Spatial  variability  in  the  near-surface  TKE  held 
has  two  distinguishable  patterns.  The  relatively  large  values  in  the  region  south  of  Cape 
Blanco,  with  a  maximum  close  to  the  shelf-break  around  42  °N,  have  a  spatial  pattern 
closely  related  to  that  of  the  wind  stress  (Figure  13)  and  clearly  represent  the  response 
of  the  surface  layer  TKE  to  the  increased  wind  stress  in  that  region.  On  the  other  hand, 
the  zonally-oriented  deep-red  patch  along  42  °N  represents  the  TKE  associated  with  the 
separated  jet.  The  origin  of  the  “jet-born”  TKE  pattern  can  be  found  in  the  area  of 
negative  vorticity  north  of  the  narrow  strip  of  zero  vorticity  along  the  jet  axis  (Figure  16, 
left).  The  negative  vorticity  appears  to  affect  the  ageostrophic  velocity,  evidently  con¬ 
tributing  to  an  increase  in  its  meridional  component  va  close  to  the  surface.  This  leads 
to  increased  vertical  shear  in  va  in  the  surface  layer,  visible  in  the  offshore  section  along 
126  °W  (Figure  14i),  and  corresponding  increased  production  of  TKE  (Figure  15g).  The 
local  increase  in  the  TKE  offshore  just  north  of  the  jet  axis  consequently  appears  to  be 


forced  in  response  to  increased  shear  in  the  cross-jet  ageostrophic  velocity.  Note  that 
the  dominant  spatial  variability  in  the  near-surface  TKE  and  hence  in  the  depth  of  the 
surface  boundary  layer  (SBL)  is  found  on  scales  smaller  than  the  scales  of  the  wind  stress. 
This  finding  suggests  that  the  approach  to  obtain  estimates  of  surface  currents  from  the 
satellite  information  alone  (SSH  and  wind  stress),  based  on  the  assumption  about  the 
spatially  uniform  SBL  depth  [e.g.,  Saraceno  et  al.,  2008],  can  be  limited  in  the  CTZ. 

To  assess  the  spatial  extent  of  unstable  stratification  indicated  by  the  vertical  sec¬ 
tions  of  N2  and  Kh  in  Figure  15  we  examine  a  horizontal  held  of  the  maximum  of  the 
stability  function  Gh  [e.g.  Wijesekera  et  al,  2003],  in  the  upper  25  m  on  1  August-0800 
(Figure  16,  right).  Negative  values  of  Gh  correspond  to  areas  of  stable  stratification,  while 
positive  values  correspond  to  areas  of  unstable  stratification  and  result  in  accompanying 
large  values  of  Kh-  We  note  that  there  are  extensive  regions  of  positive  Gh  in  the  CTZ 
surface  layer  and  that  most  of  the  patches  of  larger  horizontal  scale  are  offshore  of  126  °W. 
In  general,  the  patches  of  positive  Gh  (Figure  16,  right)  can  be  identified  with  areas  where 
the  surface  how  direction  is  not  aligned  with  the  surface  isopycnals  (Figure  13,  right)  and 
the  advection  of  heavier  water  over  light  water  occurs.  Note  that  these  regions  are  not 
necessarily  associated  with  the  most  energetic  features  of  the  surface  velocities  associated 
with  the  separated  coastal  jet. 

4.5  Time-dependent  behavior  in  the  separated  jet 

In  order  to  examine  time-dependent  behavior  in  the  separated  jet,  we  plot  in  Figure  17 
horizontal  helds  of  vertical  velocities  w  at  25  m  depth,  the  relative  vorticity  of  the  sur¬ 
face  currents,  and  the  surface  potential  density,  with  SSH  contours  superposed,  on  four 
days  with  variable  wind  conditions.  The  horizontal  region  of  the  plots  encompasses  the 
separated  jet  between  (127°W,  124. 5°W)  and  (41.7°N, 42.7°N).  The  days  and  the  cor¬ 
responding  wind  conditions  (Figure  8)  are  as  follows:  28  July-0000,  moderately  strong 
winds  relatively  early  in  a  several-day  southward  wind  event;  1  August-0800,  strong  winds 


about  6  days  into  the  same  southward  wind  event;  6  August-0000,  weak  winds  between 
southward  wind  events;  13  August-0000,  strong  winds  about  6  days  into  a  second  strong 
southward  wind  event.  Also  plotted  in  Figure  17  is  a  measure  of  the  flow  imbalance 

=  d(V|,  '  u)/d« 

/(|CI  +  <C2>1/2)’ 

[McWilliams,  1985;  compare  Capet  et  al,  2008b  (16)],  where  V/,  •  u  =  ux  +  vy  is  the 
horizontal  divergence  of  the  surface  velocity  and  (  =  vx  —  uy  is  the  relative  vorticity  of 
the  surface  currents.  The  term  <  £ 2  >V2,  where  the  brackets  represent  a  spatial  average 
over  the  horizontal  sub-region  shown  and  the  overbar  represents  a  time  average  over  the 
period  26  July-21  August,  is  added  to  avoid  singular  behavior  at  zeros  of  (.  Values  of 
|e|  <<  1  indicate  that  the  flow  is  essentially  in  balance,  typically  through  a  geostrophic 
or  a  more  general  gradient-wind  balance.  Values  of  e  =  0(1)  indicate  that  the  flow  is 
unbalanced. 

We  look  first  at  the  fields  for  1  August-0800,  which  correspond  to  the  same  time 
as  the  sections  in  Figures  14  and  15  and  the  horizontal  fields  in  Figures  13  and  16.  The 
w  held  (Figure  17e)  shows  relatively  small  scale  spatial  variations  along  the  jet  with 
a  wave  length  20-30  km.  These  disturbances  are  found  to  propagate  westward  in  the 
direction  of  the  jet  with  a  phase  speed  of  about  25  km  day-1  (Figure  18)  and  reflect  the 
presence  of  submesoscale  instabilities  in  the  surface  layer  of  the  jet.  From  the  sections 
in  Figure  14,  the  magnitudes  of  the  vertical  velocities  in  these  disturbances  can  reach 
values  of  60  m  day^1  at  the  near-coastal  section  and  40  m  day-1  at  the  offshore  section. 
The  corresponding  spatial  variations  in  the  surface  vorticity  held  (  (due  mostly  to  the 
vorticity  in  the  geostrophic  surface  currents),  in  the  surface  potential  density  held,  and 
in  the  SSH  held  are  relatively  smooth,  reflecting  the  larger  scale  jet  structure,  and  do 
not  give  strong  indications  of  variability  on  the  smaller  scales  of  the  vertical  velocity 
disturbances.  The  e  held  has  a  spatial  structure  clearly  related  to  that  of  the  w  held,  with 


magnitudes  reaching  around  0.5  in  the  regions  of  the  large  vertical  velocities.  Values  of 
e  of  that  magnitude  clearly  indicate  that  these  small-scale  disturbances  are  not  in  either 
geostrophic  or  gradient-wind  balance  and  that  they  are  unbalanced.  Note  that  evidence 
for  the  presence  of  the  instabilities  can  be  seen  in  the  1  August  horizontal  field  of  TKE 
(Figure  16)  through  the  small  scale  fluctuations  of  TKE  in  the  jet.  In  those  fluctuations, 
the  TKE  values  are  generally  reduced  in  the  regions  of  upwelling  velocities,  as  shown  in 
the  offshore  section  (Figures  14,  15). 

During  weak  wind  conditions  on  6  August,  the  (  held  (Figure  17j)  shows  strong 
along-jet  perturbations  with  a  wave-like  structure  and  an  along-jet  scale  of  about  50-70 
km.  These  perturbations,  also  present  in  the  surface  potential  density  held,  propagate 
along  the  jet  with  a  phase  speed  of  25-30  km  day-1.  In  contrast  to  the  behavior  on  1 
August,  the  w  held  (Figure  17i)  along  the  jet  exhibits  weaker  variability,  and  smaller 
magnitudes.  The  largest  magnitudes  of  w  are  associated  with  downwelling  and  are  found 
near  the  crests  of  the  along- front  instabilities  (Figure  17k),  where  the  curvature  of  the 
SSH  held  is  locally  large.  The  e  values  are  relatively  low  (Figure  171)  with  two  spots  of 
increased  values,  also  near  the  crests.  The  instabilities  present  on  this  day  appear  to  have 
a  different  dynamical  structure  than  those  found  during  strong  winds  on  1  August  and, 
from  the  e  helds,  to  be  essentially  balanced. 

During  the  second  strong  wind  event  on  13  August  the  small-scale,  unbalanced 
instabilities  are  again  present  with  substantial  fluctuations  in  the  w  held,  accompanied 
by  related  variations  in  e  (Figure  17p),  propagating  along  the  jet  (Figure  18),  similar  to 
those  on  1  August.  Fluctuations  in  the  surface  vorticity  held  and  in  the  surface  potential 
density  held,  however,  are  larger  than  those  on  1  August,  with  some  resemblance  to  those 
present  during  weak  winds  on  6  August.  Two  days  after  this  event  on  15  August,  the 
small-scale,  unbalanced  instabilities  decay  (not  shown  here)  and  the  vorticity  and  density 
holds  become  very  much  like  those  on  6  August,  with  similar  disturbances  of  wave-like 
structure.  Thus,  there  are  some  indications  that  the  small-scale,  unbalanced  instabilities 


may  subsequently  evolve  into  the  larger  scale,  balanced  instabilities. 

In  order  to  examine  the  generation  and  propagation  of  the  small-scale,  unbalanced 
instabilities  shown  in  Figure  17e,m,  we  construct  a  H ovmbller  diagram  of  the  time  de¬ 
pendent  vertical  velocities  w  at  25  m  depth,  sampled  along  the  mean  path  of  maximum 
density  gradient  \VhP\  representing  the  jet  axis,  for  two  time  intervals  when  those  insta¬ 
bilities  are  present:  28  July-9  August  and  10-22  August  (Figure  18).  During  the  first 
time  interval,  strong  instabilities,  identified  by  large  magnitude  vertical  velocities  (up  to 
60  m  day'1),  start  to  develop  on  30  July.  They  propagate  along  the  jet  path  and  stay 
in  the  subdomain  until  4  August.  Most  of  instabilities  during  that  interval  originate  be¬ 
tween  50  and  100  km  from  the  start  point  of  the  jet  path  in  a  region  where  the  jet  is  still 
mainly  flowing  toward  the  south  before  turning  offshore  (see  the  plot  of  the  mean  path  in 
Figure  18).  If  we  consider  individual  perturbations,  then  the  lifetime  of  a  perturbation, 
defined  as  the  time  when  the  magnitude  of  w  is  greater  than  25  m  day'1,  can  be  estimated 
to  be  around  3  days  and  the  distance  traveled  to  be  around  75  km,  which  gives  an  average 
propagation  speed  of  25  km  day'1.  During  the  second  time  period  of  instability  occur¬ 
rence,  more  instabilities  originate  closer  to  the  start  point  of  the  jet  path,  and  not  all  of 
them  travel  all  the  way  to  the  offshore  part  of  the  subdomain.  Instabilities  are  present  in 
the  subdomain  from  11  to  17  August.  They  are  characterized  by  larger  vertical  velocities 
up  to  100  m  day'1,  longer  lifetimes  around  4  days,  and  by  smaller  average  propagation 
speeds  of  18  km  day'1. 

To  further  examine  the  time-dependent  relationship  of  the  relatively  large  fluctua¬ 
tions  in  vertical  velocity,  that  contribute  to  the  horizontal  fields  of  mean  and  standard 
deviation  of  w  in  Figure  10,  with  the  occurrence  of  the  submesoscale  instabilities,  we 
calculate  the  spatial  average  (again  denoted  by  brackets)  of  the  vertical  velocity  squared 
<  w2  >  at  25  m  depth  and  surface  vorticity  squared  <  (2  >  over  the  sub-region  around 
the  offshore  separated  jet  shown  in  Figure  17.  Those  values,  together  with  the  magni¬ 
tude  of  wind  stress,  are  plotted  as  time  series  for  26  July-22  August  in  Figure  19.  Two 


events  of  large  <  w2  >  occur,  on  3  and  13  August,  respectively.  These  events  coincide 
with  the  presence  of  the  small-scale,  unbalanced  instabilities.  The  lagged  relation  of  the 
larger  values  of  <  w2  >,  with  larger  values  of  the  southward  wind  stress  is  clear.  (The 
maximum  correlation  coefficient  of  <  w2  >  and  \tv\  is  0.62  at  a  lag  of  2  days.)  The  two 
events  of  large  <  w2  >  are  followed  by  increased  <  (2  >  with  about  2  days  lag.  This 
evidently  reflects  the  weakening  of  the  small-scale,  unbalanced  instabilities  and  the  sub¬ 
sequent  growth  of,  or  evolution  into,  the  larger  scale  balanced  instabilities,  characterized 
by  substantial  fluctuations  in  vorticity. 

4.6  Frontogenesis  along  the  separated  jet 

In  order  to  assess  the  frontogenesis  along  the  separated  jet  we  use  frontogenesis  function 
F  following  Capet  et  al,  2008b: 


D\Vhp\2 

Dt 


2  F, 


(6) 


where  |V/,p|2  =  (d p/dx)2  +  ( dp/dy )2  is  the  absolute  value  of  surface  density  gradient 
squared.  Horizontal  fields  of  the  kinetic  energy  of  surface  currents  KE  =  0.5(w2  +  v2), 
iV/jpl2,  and  the  part  of  F  due  to  horizontal  advection  [ Hoskins ,  1982;  Capet  et  al.,  2008b]: 
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averaged  in  time  over  the  27  day  period  26  July-21  August  are  shown  in  Figure  20.  The 
structure  of  the  surface  kinetic  energy  held  is  remarkable  in  that  by  far  the  largest  values 
over  the  entire  domain  are  concentrated  in  the  region  of  the  separated  jet  offshore  south  of 
Cape  Blanco.  Dynamically  consistent  large  values  of  the  absolute  value  of  the  horizontal 
gradient  of  the  surface  density  squared  |V/;p|2  are  also  found  concentrated  along  the  jet 
in  the  same  region.  The  largest  values  of  \Vhp\2  occur,  however,  near  the  coast  between 


42.5  °N  and  43  °N,  where  the  jet  is  still  flowing  south  in  the  general  direction  of  the 
wind  stress.  The  large  values  of  |  V *yo|2  in  the  separated  jet  are  clearly  related  to  the 
concentrated  large  mean  values  of  the  frontogenesis  function  Fs  found  in  the  same  region. 
In  particular,  there  is  an  obvious  match  between  the  increased  magnitudes  of  both  fields  in 
the  jet  near  the  coast  between  42.5  °N  and  43  °N.  It  appears  that  frontogenesis  processes 
related  to  the  CTZ  flow  fields  play  a  major  role  in  establishing  the  structure  and  intensity 
of  the  separated  jet.  In  the  near-shore  southward  flowing  part  of  the  jet,  frontogenesis 
processes,  intensified  by  the  along-jet  wind  stress,  contribute  directly  to  the  development 
of  the  FSC.  In  the  offshore  part  of  the  jet,  frontogenesis  processes  influence  the  structure 
of  the  jet  and  thus  affect  the  behavior  of  the  submesoscale  instabilities.  Consequently, 
in  both  parts  of  the  jet,  frontogenesis  processes  play  a  role  in  the  maintenance  of  the 
concentrated  large  mean  and  RMS  vertical  velocities  shown  in  Figure  10. 

To  document  one  additional  aspect  of  the  behavior  of  the  vertical  velocities,  regard¬ 
ing  possible  differences  in  the  maximum  magnitudes  of  the  downwelling  and  upwelling 
velocities,  we  calculate  probability  density  functions  (PDFs)  of  the  extreme  positive  and 
negative  values  of  w  over  the  top  80  m  depth  and  over  the  time  period  26  July-21  August 
(Figure  21).  Two  horizontal  regions  are  selected,  one  for  the  near-coastal  part  of  the  jet, 
where  it  is  flowing  southward  in  the  direction  of  the  wind  stress,  and  one  for  the  offshore 
part,  where  it  flows  westward.  The  near-coastal  region  extends  from  42.3  °N  —  43  °N  and 
includes  all  values  along  east-west  grid  lines  18  km  on  either  side  of  the  jet  axis,  defined 
by  the  mean  location  of  the  maximum  surface  velocity.  The  offshore  section  extends 
from  125.7 °W  —  126.3  °N  and  includes  all  values  along  north-south  grid  lines  18  km  on 
either  side  of  the  jet  axis.  From  Figure  21a  we  can  see  that  in  the  near-coastal  region, 
where  the  down-front  winds  evidently  contribute  to  an  increase  in  the  absolute  values  of 
surface  density  gradients  (Figure  20)  and  to  an  intensification  of  the  accompanying  FSC 
(Figure  10),  the  positive  upwelling  velocities  clearly  have  more  occurrences  of  larger  mag¬ 
nitude  than  do  the  negative  downwelling  velocities.  Thus,  in  this  region  we  find  typically 


stronger  upwelling  velocities  and  weaker  downwelling  velocities.  Consistent  with  that  re¬ 
sult,  stronger  upwelling  velocities  are  evident  in  the  instantaneous  near-coastal  section 
in  Figure  14.  The  behavior  differs,  however,  in  the  offshore  region  (Figure  21b),  where 
the  downwelling  velocities  have  more  occurrences  of  larger  magnitude  than  the  upwelling 
velocities.  Thus,  in  the  offshore  region,  where  the  wind  stress  is  normal  to  the  direction 
of  the  jet  and  where  the  more  energetic  vertical  circulation  appears  to  be  primarily  as¬ 
sociated  with  time-dependent,  submesoscale  instabilities,  we  find  evidence  for  typically 
stronger  downwelling  velocities.  That  general  behavior  is  also  consistent  with  the  model 
results  of  Capet,  et  al.  [2008b]  and  with  theoretical  results  [ Hoskins ,  1982]  concerning  the 
FSC  associated  with  frontogenesis  caused  by  larger  scale  horizontal  deformation  fields. 

5  Discussion 

The  phenomenon  of  nonlinear  interaction  between  effects  of  a  geostrophic  flow  and  a  wind 
stress  applied  to  a  surface  Ekman  layer  has  been  studied  in  a  number  of  papers  [Stern, 
1965;  Niiler,  1969;  Thomas  and  Rhines,  2002;  Thomas  and  Lee ,  2005;  Pedlosky,  2008].  Lee 
et  al.  [1994]  used  a  two-dimensional  numerical  model  to  study  the  case  when  a  spatially 
uniform  wind  is  applied  both  along  and  across  a  geostrophically  balanced  jet.  They 
reported  secondary  circulation  that  is  created  as  a  result  of  a  nonlinear  interaction  between 
the  jet  and  the  wind  driven  flow  in  the  Ekman  layer.  In  that  study  (as  well  as  studies 
mentioned  above)  an  expression  similar  to  equation  (4)  for  the  Ekman  transport  Me 
was  used.  They  found  that  when  the  wind  blows  perpendicular  to  the  jet  the  secondary 
circulation  is  50  %  weaker  than  that  when  the  wind  blows  parallel  to  the  jet,  with  upward 
vertical  advection  on  the  upwind  side  of  the  jet  (C  <  0)  and  downward  advection  on 
downwind  side  (£  >  0).  Lee  and  Niiler  [1998]  modeled  the  ocean  response  to  uniform  wind 
stress  forcing  over  geostrophically  balanced  eddies.  Their  results  for  the  mean  secondary 
circulation,  were  consistent  with  the  features  found  in  Lee  et  al.  [1994],  Centurioni  et 
al.  [2008]  used  surface  drifter  data,  satellite  SSH  measurements,  and  NCEP  reanalysis 


winds  to  map  the  time-average  15  m  depth  geostrophic  velocity  field  in  the  CCS.  The 
resultant  mean  circulation  and  eddy  energy  distributions  were  found  to  be  in  reasonable 
agreement  with  ROMS  CCS  model  results  from  Marchesiello  et  al.  [2003].  Analysis 
of  the  time-average  near-surface  ageostrophic  velocity  field  in  the  ROMS  CCS  solutions 
showed  behavior  similar  to  that  found  in  Lee  and  Niiler  [1998]  and  thus  was  argued 
to  be  associated  with  the  nonlinear  interaction  of  Ekman  dynamics  with  the  geostrophic 
vorticity  field.  Pedlosky  [2008],  in  an  analytical  study  for  a  homogeneous  fluid,  determined 
the  nonlinear  effects  of  geostrophic  vorticity  ( g  and  wind  stress  curl  on  the  Ekman  layer 
thickness.  The  importance  of  accounting  for  the  Ekman  layer  depth  for  calculations  of 
surface  velocity  from  satellite  SSH  measurements  was  noted  in  Section  4.4.  Thomas  and 
Lee  [2005]  studied  the  effect  of  strong  winds  blowing  in  the  direction  of  a  frontal  jet 
on  frontogenesis  and  found  that  the  frontogenesis  secondary  circulation,  characterized 
by  subduction  on  the  dense  side  of  the  front  and  upwclling  along  the  frontal  interface,  is 
intensified.  The  results  of  that  study  appear  especially  relevant  to  the  behavior  found  here 
in  the  near-coastal  region  of  the  separated  jet  where  it  flows  southward  in  the  direction  of 
the  wind  stress.  Our  simulations  of  the  3D  dynamics  of  the  jet  in  Oregon  CTZ,  however, 
contain  a  more  complicated  set  of  time-  and  space-dependent  dynamical  processes  than 
those  represented  in  most  of  these  idealized  model  investigations.  Additional  process- 
oriented  studies  in  a  realistic  CTZ  environment  would  be  useful. 

The  behavior  and  the  dynamics  of  circulation  processes  near  the  ocean  surface  have 
been  addressed  further  in  several  recent  studies  [Capet  et  al .,  2008a, b,c;  Lapeyre  and  Klein, 
2006;  Lapeyre  et  al.,  2006;  Klein  et  al.,  2008;  Boccaletti  et  al.,  2007;  Fox-Kemper  et  al., 
2008;  Fox-Kemper  and  Ferrari,  2008].  In  particular,  the  subject  of  near-surface  vertical 
exchange  associated  with  submesoscale  processes  has  been  discussed  in  a  useful  recent 
review  paper  by  Klein  and  Lapeyre  [2009].  The  previous  results  of  most  relevance  here 
are  those  of  Capet  et  al.  [2008a, b,c]  where  the  dynamics  of  near  surface  submesoscale  pro¬ 
cesses  in  an  idealized  eastern  boundary  current  oceanic  regime  with  steady  wind  forcing 


was  investigated.  In  particular,  Capet  et  al.  [2008b]  found  evidence  for  more  intense  FSC 
in  fronts  oriented  in  the  downwind  direction,  consistent  with  the  results  of  Thomas  and 
Lee  [2005]  and  in  agreement  with  the  behavior  found  here  in  the  near-coastal  southward 
flowing  part  of  the  separated  coastal  jet.  Capet  et  al.  [2008b]  also  found  evidence  for 
frontogenesis  processes  associated  with  regions  of  large  surface  density  gradients,  similar 
to  the  behavior  shown  here  in  Figure  20.  In  addition,  Capet  et  al.  [2008b]  found  occur¬ 
rences  of  submesoscale  instability  growth  and  propagation  along  fronts  with  instability 
wavelengths  of  20-30  km.  The  instabilities  were  characterized  by  spatial  variations  in  the 
surface  potential  density  erg  and  by  increases  in  magnitude  of  corresponding  submesoscale 
surface  vorticity  Q  fluctuations.  The  small-scale,  unbalanced  instabilities  found  here  (Fig¬ 
ures  17,18)  have  similar  spatial  scales,  but  appear  to  be  characterized  most  strongly  by 
substantial  fluctuations  in  w,  resulting  in  vigorous  vertical  exchange  in  the  surface  layer 
(Figure  14j),  with  relatively  smaller  variations  in  (.  Moreover,  there  is  a  clear  relation 
between  increased  strength  of  the  wind  stress  and  occurrence  of  these  instabilities  (Fig¬ 
ure  19).  To  what  extent  the  development  of  these  instabilities  depend  on  the  wind  stress 
through  resultant  changes  in  the  structure  of  the  separated  jet,  or  on  the  increase  in 
magnitude  of  the  surface  layer  ageostrophic  transport,  or  on  their  interaction,  is  not  clear 
from  these  simulations.  Additional  controlled  process  studies  to  address  these  questions 
would  be  useful.  In  contrast,  the  somewhat  larger  scale  instabilities  found  on  6  August 
are  characterized  by  substantial  variations  in  surface  (  and  erg  and  by  a  larger  degree  of 
balance  (Figure  17j,k,l). 

We  note  that  the  primary  differences  in  the  results  of  our  limited-time  simulation  of 
late  summer  conditions,  with  realistic  wind  stress  forcing  and  realistic  coastal  topography, 
and  those  from  the  longer-time  simulations  of  Capet  et  al.  [2008a, b,c],  with  more  idealized 
forcing  and  topography,  is  in  the  strong  space-  and  time-dependent  behavior  found  here. 
In  particular,  the  significant  spatial  concentration  of  the  near-surface  ageostrophic  vertical 
circulation  processes  in  the  region  of  the  separated  coastal  jet  off  Cape  Blanco  (Figure  10) 


and  the  time  dependent  behavior  of  that  circulation  in  response  to  the  wind  stress  forcing 
(Figure  19),  appears  to  identify  potentially  important  components  of  realistic  CTZ  flow 
fields  off  Oregon  and  northern  California. 

6  Summary 

Based  on  comparisons  with  different  types  of  data,  our  3D  nonlinear  model  describes  the 
dynamics  both  on  the  shelf  and  in  the  CTZ  off  Oregon  qualitatively  correctly.  Com¬ 
parisons  with  currents  at  mid-shelf  moorings  indicate  that  the  model  reproduces  both 
depth- averaged  and  baroclinic  dynamics  on  the  shelf.  The  modeled  SST  front  evolution 
is  consistent  with  satellite  SST  fields  showing  the  major  feature  of  summer-time  Oregon 
CTZ  dynamics,  such  as  the  separation  of  the  coastal  jet  near  Cape  Blanco.  HF  radar 
measurements  of  surface  currents  (Figure  6)  provide  evidence  that  the  observed  jet  that 
separates  off  Cape  Blanco  is  as  energetic  as  the  modeled  separating  jet  with  maximum 
speeds  reaching  0.6-0. 8  m  s-1.  Lagrangian  analysis  reveals  characteristic  flow  patterns 
over  the  shelf  and  in  the  CTZ  off  Oregon,  in  particular,  separated  jets  near  Heceta  Bank 
and  Cape  Blanco. 

The  behavior  of  near-surface  circulation  in  the  CTZ  off  Oregon  during  a  late-summer 
27  day  time  period,  26  July-21  August,  was  investigated.  Significant  vertical  velocities 
were  found  concentrated  in  the  region  of  the  separated  jet  south  of  Cape  Blanco.  In  the 
near-coastal  part  of  the  separated  jet,  where  it  was  flowing  southward  in  the  direction  of 
the  winds,  the  wind  stress  appeared  to  intensify  the  frontogenesis  secondary  circulation. 
The  associated  instantaneous  vertical  velocities  were  as  large  as  100  m  day-1.  In  the 
offshore  part  of  the  separated  jet,  where  it  was  flowing  westward,  the  energetic  vertical 
velocities  were  primarily  associated  with  time-dependent  submesoscale  instabilities.  Those 
instabilities  were  present  during  strong  southward  wind  events.  They  propagated  in  the 
direction  of  the  jet  with  wavelengths  20-30  km  and  propagation  speeds  of  18-25  km  day-1. 
The  larger  vertical  velocities  were  typically  50  m  day-1. 


As  a  next  step,  it  would  be  useful  to  extend  these  analyses  to  a  coupled  bio-physical 
model,  to  see  how  the  surface  layer  processes  within  the  jet  affect  biological  variability. 
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Figure  1:  Monthly  SST  composites  from  GOES  satellite  (top)  and  monthly  averaged  SST 
from  ROMS  (bottom)  for  May-August  (left  to  right).  Coastal  HF  radar  locations  are 
shown  in  the  lower  right  panel  by  blue  circles,  and  mid-shelf  moorings  by  magenta  circles. 
’OB’  denotes  Cape  Blanco.  The  black  line  shows  the  200  m  isobath. 


Figure  2:  Along-shore  COAMPS  wind  stress  Ty  at  the  NH10  (44.6  °N)  (a)  and  Rogue  River 
(42.4  °N)  (b)  moorings;  and  depth-averaged  along-shore  current  V  at  NH10  (44.6  °N)  (c), 
Coos  Bay  (43.2  °N)  (d),  and  Rogue  River  (42.4  °N)  (e)  mooring  observations  (gray  lines) 
and  model  (black  lines).  Variables  are  presented  along  respective  major  principal  axes. 
Root  mean  square  error  (RMSE)  and  correlation  coefficients  (CC)  for  data  and  model  are 
shown. 


Figure  3:  Mean  and  variance  ellipses  for  depth-averaged  currents  at  the  mooring  locations 
(top  to  bottom:  NH10,  Coos  Bay,  Rogue  River)  for  observations  (gray),  and  model  (black) 
over  the  calculation  interval  (see  Figure  2). 


Figure  4:  Vertical  profiles  of  time-averaged  statistics  for  cross-shore  u  (pale  colors, 
thick  lines)  and  along-shore  v  (bright  colors,  thin  lines)  components  of  velocity  pro¬ 
jected  on  respective  principal  axes  of  the  depth- averaged  velocities:  means  (a),  stan¬ 
dard  deviations  (std)  (b),  normalized  root  mean  square  error  (NRMSE)  (c),  defined  as 

1 

NRMSE  =  (u0bs  —  Umod)2 / u2obs  ,  where  the  over-bar  denotes  a  time-average,  corre¬ 
lation  coefficients  (CC)  (d).  In  (a)  and  (b),  red  lines  denote  observations,  blue  lines  - 
model  results.  Data  is  taken  from  mid-shelf  moorings:  NH10  (top),  Coos  Bay  (middle), 
and  Rogue  River  (bottom). 


Figure  5:  Surface  current  statistics  from  HF  radar  and  model  for  May-August  (top  to 

bottom):  two  left  panels  -  monthly  averaged  current  vectors  u  and  speed  |u|  in  color; 

1 

-  2 

two  right  panels  -  RMS  speed  deviations  from  the  mean  Urms  —  \u  ~  h)2  +  (v  —  v)2  , 

where  the  over-bar  denotes  a  time-average.  Color  contour  intervals  for  Urms  are  0-03 
m  s^1.  The  white  contour  on  the  model  fields  shows  the  area  of  the  corresponding  data 
coverage. 

Figure  6:  Instantaneous  fields  of  surface  current  vectors  u  and  speed  |u|  (color)  for  HF 
radar  (top  panel)  and  model  (middle  panel)  in  July-August  in  the  area  near  Cape  Blanco, 
and  observed  satellite  SSH  (from  track  206)  (red)  and  model  SSH  sampled  at  the  same 
location  (blue)  shown  in  the  bottom  panels.  The  SSH  has  the  mean  taken  out.  The 
location  of  track  206  is  denoted  by  a  white  dashed  line. 

Figure  7:  Cross-shore  sections  of  potential  density  erg  measured  during  a  SeaSoar  survey 
(a,c,e),  and  model  fields  (b,d,f)  sampled  at  the  same  times  and  locations  as  the  observa¬ 
tions,  along  44.25  °N  during  2  August  (a,b),  43.5  °N  during  3-4  August  (c,d)  and  41.9  °N 
during  6-7  August  (e,f).  Color  contour  intervals  are  0.25  kg  m”3. 

Figure  8:  Time  series  of  the  north-south  component  of  the  wind  stress  Ty,  averaged  over 
the  across-jet  section  along  126  °W  between  41.65  °N  and  42.4  °N  (see  Figure  13),  from  26 
July-31  August  2002.  Black-shaded  bars  denote  strong  wind  events,  gray-shaded  bars  - 
weak  wind  events. 

Figure  9:  Surface  potential  density  cre  in  color  and  SSH  in  white  contours  (left)  and 
surface  current  vectors  u  and  their  standard  deviations  in  color  (right),  averaged  over  26 
July-21  August.  Contour  intervals  for  SSH  are  0.03  m.  The  solid  black  line  shows  the 
200  m  isobath. 

Figure  10:  Mean  vertical  velocity  w  at  25  m  depth  (left)  and  its  standard  deviation  (right) 
in  color  and  SSH  in  black  contours,  averaged  over  26  July-21  August.  Contour  intervals 
for  SSH  are  0.03  m.  The  solid  black  line  shows  the  200  m  isobath. 


Figure  11:  Lagrangian  surface  particle  trajectories  (gray  lines)  during  1-10  August  (left) 
and  1-20  August  (right).  Particles  are  released  simultaneously  on  1  August  along  the 
81  m  isobath  (depth  of  NH10  mooring)  every  10  km  in  the  meridional  direction.  Circles 
denote  particles  final  locations  if  they  are  within  the  domain.  The  coastline  is  shown  in 
bold,  isobaths  100,  200,  and  1000  m  -  in  thin  black  lines. 


Figure  12:  Divergence  (a,c)  and  deformation  (b,d)  terms  of  the  surface  strain  rate  field 
(3)  averaged  over  26  July-21  August  (a,b)  and  on  1  August-0800  (c,d).  The  solid  black 
line  shows  the  200  m  isobath. 


Figure  13:  Surface  wind  stress  magnitudes  in  color  and  directions  in  vectors  (left),  surface 
potential  density  in  color  and  surface  velocity  u  in  vectors  (right)  on  1  August-0800.  Cross¬ 
jet  section  locations  along  42.63  °N  and  126  °W,  respectively,  are  shown  (right)  by  straight 
lines  where  the  middle  tick  denotes  the  jet  axis  location.  The  solid  black  line  shows  the 
200  m  isobath. 


Figure  14:  Cross-jet  near-coastal  (a-e)  and  offshore  (f-j)  sections  on  1  August-0800:  po¬ 
tential  density  (a,f),  geostrop  hie  vg  (b),  ug  (g)  and  ageostrophic  va  (c),  ua  (h)  components 
of  along-jet  velocity,  ageostrophic  component  of  the  cross-jet  velocity  ua  (d),  va  (i)  and 
vertical  velocity  w  (e,j).  Color  contour  intervals  are  0.1  kg  m-3  for  erg,  0.05  nr  s_1  for 
(■ ug,vg ),  0.025  m  s'1  for  ( ua,va )  and  5  m  day-1  for  w.  Bold  isolines  denote  zero  veloc¬ 
ity.  The  near-coastal  cross-jet  section  is  along  42.63  °N,  and  the  offshore  section  is  along 
126 °W  (see  Figure  13). 


Figure  15:  Cross-jet  near-coastal  (a-d)  and  offshore  (e-h)  sections  on  1  August-0800: 
buoyancy  frequency  N2  (a,e),  vertical  turbulent  diffusivity  coefficient  Kh  (b,f),  turbulent 
kinetic  energy  in  color,  Richardson  number  Ri  =  N2 / S2  in  white  contours  (c,g)  and  shear 
production  P  =  KmS2  (d,h).  Negative  values  of  N2  are  shown  with  white  stars.  The 
near-coastal  cross-jet  section  is  along  42.63  °N,  and  the  offshore  section  is  along  126  °W 
(see  Figure  13). 

Figure  16:  Ageostrophic  component  of  the  surface  current  ua  in  vectors  and  relative 
vorticity  of  geostrophic  surface  current  ( g  in  color  (left),  maximum  of  turbulent  kinetic 
energy  in  upper  25  m  in  color  and  SSH  in  contours  (middle),  maximum  of  the  stability 
function  (S'/,  in  upper  25  m  (right)  on  1  August-0800.  The  solid  black  and  white  lines 
show  the  200  m  isobath. 

Figure  17:  Vertical  velocity  at  25  m  depth  w  (a,e,i,m),  relative  vorticity  of  the  surface 
current  (  (b,f,j,n),  surface  potential  density  erg  (c,g,k,o)  and  a  measure  of  the  imbalance 
e  (d,h,l,p)  (see  (5)  Section  4.5)  for  the  area  around  the  separated  jet  on  28  July-0000 
(a-d),  1  August-0800  (e-h),  6  August-0000  (i— 1) ,  13  August-0000  (m-p).  Black  and  white 
contours  show  the  SSH  with  contour  intervals  0.02  m. 

Figure  18:  Hovmbller  diagrams  of  w  at  25  m  sampled  along  the  jet  path  (shown  in  small 
windows),  as  determined  from  the  mean  location  of  the  maximum  in  |V/,p|,  for  two  time 
intervals  of  instability  occurrence. 

Figure  19:  Time  series  plots  of  the  spatial  average,  denoted  by  brackets,  over  the  region 
around  the  separated  jet  shown  in  Figure  17  of  <  w2  >  at  25  m  (m2  s~2xl0~8)  (red  line) 
and  surface  <  (2  >  (s_2x5xl0^10)  (blue  line).  Also,  plotted  is  the  corresponding  time 
series  of  \rv\  (N  m-2)  (gray  line)  calculated  as  in  Figure  8. 

Figure  20:  Kinetic  energy  of  surface  current  (left),  absolute  value  of  the  horizontal  gradient 
of  surface  density  squared  | V*,p|2  (middle)  and  the  frontogenesis  function  Fs  (7)  (right), 
averaged  over  26  July-21  August.  The  solid  black  line  shows  the  200  m  isobath. 


Figure  21:  Histogram  of  the  extreme  positive  and  negative  values  of  the  vertical  velocity 
w  over  the  top  80  m  depth  across  the  coastal  jet  in  its  near-coastal  (a)  and  offshore  (b) 
parts.  The  near-coastal  part  of  the  jet  for  this  sampling  extends  from  42.2  °N  —  43  °N  and 
the  offshore  part  extends  from  125.7  °W  —  126.3  °W.  The  sampling  is  made  in  a  36-km 
band,  centered  along  the  jet  axis,  defined  here  by  the  mean  location  of  the  maximum 
surface  velocity,  and  the  integration  interval  is  from  26  July-21  August. 
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